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Abstract 



We have considered a cosmological model of the FRW universe with variable 
G and A. The solutions have been obtained for flat model with particular form 
of cosmological constant. The cosmological parameters have also been obtained 
for dust, radiation and stiff matter. The statefinder parameters are analyzed and 
have shown that these depends only on w and e. Further the lookback time, 
proper distance, luminosity distance and angular diameter distance have also been 
calculated for our model. 
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I. INTRODUCTION 

The Einstein field equation has two parameters - the gravitational constant G and the 
cosmological constant A. The Newtonian constant of gravitation G plays the role of a 
coupling constant between geometry and matter in the Einstein field equations. In an 
evolving Universe, it appears natural to look at this "constant" as a function of time. 
Numerous suggestions based on different arguments have been proposed in the past few 
decades in which G varies with time [H [2]. Dirac [3] proposed a theory with variable G 
motivated by the occurrence of large numbers discovered by Weyl, Eddington and Dirac 
himself. Many other extensions of Einstein's theory with time dependent G have also been 
proposed in order to achieve a possible unification of gravitation and elementary particle 
physics or to incorporate Mach's principle in general relativity [I]. From the point of view 
of incorporating particle physics into Einstein's theory of gravitation, the simplest approach 
is to interpret the cosmological constant A in terms of quantum mechanics and the Physics 
of vacuum [3]. The A term has also been interpreted in terms of the Higgs scalar field [BJ. 

A as a function of time has also been considered by several authors in various variable 
G theories in different contexts [7] . With this in view, several authors [H [9] have proposed 
linking the variation of G with that of A in the framework of general relativity. By 
considering the conservation of the energy-momentum tensor of matter and vacuum take 
together, many authors have invoked the idea of a decreasing vacuum energy and hence 
a varying cosmological constant A with cosmic expansion in the frame work of Einstein's 
theory. Present-day astronomical observations indicate that the cosmological constant 
A is extremely negligible being < 10~ 56 cm 2 . But the value of this constant should be 
10 50 times larger according to the Glashow-Weinberg-Salam model [11] for electro-weak 
unification and 10 107 times larger according to GUT [12] for grand unification. 

In attempt to modify the general theory of relativity, Al-Rawaf and Taha [13] related 
the cosmological constant to the Ricci scalar R. This is written as a built-in-cosmological 
constant, i.e., A oc R. Since the Ricci scalar contains a term of the form -, one adopts this 
variation for A i.e., A oc - [14] . Similarly, another two forms for A have been chosen as: 
A oc p and A oc ^ [UJ; where p is the energy density. 
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II. THE MODEL 



The action of our model is 

S I d 4 xC 



(1) 



where G is the Newton's gravitational constant and F(G) is an arbitrary function of G while 
C m is a matter Lagrangian. From the Euler-Lagrange equation 



dC „ dC 

— = V, 
dG 



d(d,G) 



to obtain 



<9F _ 
~dG ~ CP' 



Moreover from the variation with respect to g^ u , we have from ([I]), 



Writing 



we get from Q 



-GF(G) = A(t), 



R^u - \g^R = 8nGT^ u + g^„A(t). 



(2) 



(3) 



(4) 



(5) 



(6) 



As a background geometry, we consider a spatially homogeneous and isotropic FRW line 
element 



d s 2 = dt 2 - a 2 (t) 



dr 



+ r 2 (d6 2 + sin 2 



(7) 



.1 — kr 2 

where a(t) is the scale factor and k = — 1,0, +1 is the curvature parameter for spatially 
open, flat and closed universes, respectively. 

We assume that the cosmic matter is represented by the energy-momentum tensor of a 
perfect fluid 



T, 



fJ,V 



{p + p)u^u u -pg^. 



(8) 



Here p and p are the energy density and pressure of the cosmic matter, is the four velocity 
satisfying = 1. 
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We take the equation of state 



P 



{w-l)p, l<w<2. 



(9) 



Now the for the metric (J71) and stress-energy tensor ([8]), the field equations (|6j) take the form 

2 k 



-2- 



+ A = 87rGp, 



3^V + 3 A_A = 8 7rG'p. 
va/ a 2 

In view of vanishing of divergence of Einstein tensor, we have 



87rG(p + 3(p + p)-J + 8irpG + A = 0. 
The usual energy conservation equation = yields 



p + 3(p + p)- = 0. 
a 



Using (13) in (12) gives 



8npG + A = 0. 



(10) 
(11) 

(12) 

(13) 
(14) 



Note that (14) implies that G is constant whenever A is constant and vice- versa. Making 



use of (\9h in (14), we obtain 



= Ci_ 

P a 3w- 



(15) 



where C\ is a constant of integration. We can determine the value of C\ by assuming w = wq, 
p = p c due to spatial flatness (p c being the critical density) at t — to, where subscript '0' 



corresponds to the present value. It turns out G\ = p c a^° . Thus (|15j) becomes 



,3u;o 



p = p. 



c a 3w 



(16) 



For spatially flat case k = 0, (10) and (11) become 



8nGp = H 2 (2q- 1) + A, (17) 

87rGp = 3H 2 - A, (18) 

where H = - is the Hubble parameter, q = — 1 — is the deceleration parameter and 
G = 3H is the expansion scalar. 

Sahni et al [T7j introduced a pair of cosmological diagnostic pair {r, s} which they termed 
as Statefinder. The two parameters are dimensionless and are geometrical since they are 
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derived from the cosmic scale factor alone, though one can rewrite them in terms of the 
parameters of dark energy and matter. Additionally, the pair gives information about dark 
energy in a model independent way i.e. it categorizes dark energy in the context of back- 
ground geometry only which is not dependent on the theory of gravity. Hence geometrical 
variables are universal. Also this pair generalizes the well-known geometrical parameters 
like the Hubble parameter and the deceleration parameter. This pair is algebraically related 
to the equation of state of dark energy and its first time derivative. 

The statefinder parameters were introduced to characterize primarily flat universe (k = 0) 
models with cold dark matter (dust) and dark energy. They were defined as 

r ■ SJ- (19 > 

T — 1 

* = q7 XT- (2°) 

3(9- 5) 

For cosmological constant with a fixed equation of state (w = —1) and a fixed New- 
ton's gravitational constant, we have {1,0}. Moreover {1,1} represents the standard cold 
dark matter model containing no radiation while Einstein static universe corresponds to 
{00, —00} [IS] . In literature, the diagnostic pair is analyzed for various dark energy candi- 
dates including holographic dark energy [Tj5] , agegraphic dark energy [20] , quintessence [21] , 
dilaton dark energy [22], Yang-Mills dark energy [23J, viscous dark energy [24J, interacting 
dark energy [25], tachyon [26], modified Chaplygin gas [27] and f(R) gravity [28] to name a 
few. 

III. SOLUTIONS OF DYNAMICAL EQUATIONS 

We here assume the ansatz [29J for the variable cosmological constant 

A(t) = eH\ (21) 
where e is a constant. Making use of ({9]), (17), (18) and (21 ), we obtain a differential equation 



Solving (22), we get 



2H + (3 - e)wH 2 = 0. (22) 



= (.-3H + 2CV < 23 > 
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where C2 is an integration of constant. Note that C2 can be determined like before to get 

2H 



Hit) 



Solving for t we get 



H (3 - e)(wt - w t ) + 2' 
w(e — 3)H \ H J w 



From (18) and (21), we get 



P 



{3-e)H 2 . 



Integrating (24), we get 



a(t) = C 3 [(3 - e)(tw - t w )H + 2]^ 



Using (27) in (16), we get 



p(t) = Pc al W0 Cf w [{3 - e){tw - t w )H + 2} 



(£-3) 



Using (27) in (21), we get 



A(t) 



[H (3-e){tw-t w ) + 2Y 



Making use of (|29j) and (|28j) in (|14j) to get 

Hl(3 - e)[(3 - e)(tw - t w )H + 2]^ 



G(t) 

The deceleration parameter becomes 



2np c a 3 W0 C^ w 



q = -1 + 



w(3 - e) 



(24) 
(25) 
(26) 

(27) 

(28) 

(29) 

(30) 
(31) 



For an expanding universe, q < —1 which constrains e < 3. The expansion scalar takes the 
form 



e 



QH 



H (3-e)(wt- w t ) + 2' 
For flat universe, the density parameter can be obtained as 



(32) 



3H 2 3 

We now calculate the statefinder parameters {r, s} as well 

a r — 1 

r = -^7;, s 



aH 3 ' 



The parameters take the form 



3(?-|) 

)rm 



r = -[1 + w(e - 3)] [2 + u>( e - 3)], s = -(3 - e). 



Note that for e = 3, {r, s}={1,0} representing a static cosmological constant. 



(33) 



(34) 



(35) 



7 



A. Universe containing only dust 

The matter (baryonic and non-baryonic) satisfies the EoS parameter w — 1. Thus cos- 
mological parameters take the form 



a(t) 

Pit) 
A(t) 

G(t) 

1 




C 3 [(3-e)(t-t w )H + 2}^, 
p c a 3 W0 C^[(3 - e)(t - t w )H + 2] A 



AeHl 



[H (3-e)(t-t Wo) + 2}^ 

H 2 {3 - e)[(3 - e)(t - t w )H + 2]^ 



2np c a 3 WQ C 3 3 



1-e 

~2~' 



6H 



H {3-e)(t-w t ) + 2' 
l(e-2)(e-l), 



(3-e). 



(36) 
(37) 
(38) 

(39) 

(40) 
(41) 

(42) 
(43) 



B. Universe containing only radiation 



The radiation satisfies the EoS parameter w — 4/3. Thus cosmological parameters take 
the form 



ait) = c, (3 - .;/ /„«•,))//„ • 2 



2(3-6) 



p(t) = Pc al W0 C^ (3 - e) ( ^ - t ^o ) #o + 2 



(e-3) 



A(t) 



G(f) = 



4etf 2 



(3-e)(ft- W)#o + 2 

# 2 (3-e)[(3-e)(fi-* ™o)Ho + 2 
27ip c a 3 W0 C^ 



2e 
3-e 



(44) 
(45) 
(46) 

(47) 
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q = 1 - f , (48) 

e = T7 ^- , (49) 

(3-e)(jt-to^o)#o + 2 

r = l(4e-9)(2e-3), (50) 

s=^(3-c). (51) 



C. Universe containing only stiff matter 

The stiff fluid satisfies the EoS parameter w = 2. Thus cosmological parameters take the 
form 

a{t) = C 3 [(3 - e)(2t - t w )H + 2]^, (52) 

p(t) = p c a 3 w °C^[(3 - e)(2t - t w )H + 2] A, (53) 
AeH 2 

A(t) = [H (3-e)(2t-t w ) + 2]^ (54) 
G(t) = H^3-e)[(3-e)(2t-t w )H + 2]^ 

1 ' 2np c a 3 Q m C^ ' 1 ' 

q = 2 - e, (56) 

Q _ 6_^o / 57 x 

"~ H (3 - e)(2t - w t ) + 2' 1 J 

r = (2e-5)(e-2), (58) 

s = ^(3-e). (59) 



IV. SOME CONSEQUENCES 

In this section, we'll discuss lookback time, proper distance, luminosity distance and 
angular diameter as the following [30J. 



A. Lookback time 

If a photon emitted by a source at the instant t and received at the time to then the 
photon travel time or the lookback time t — to is defined by 

f a da , . 

t-t = — 60 
Ja a 
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where ao is the present value of the scale factor of the universe and can be obtained from 



(27) as (at t = t , w = w ) 

a = 2^=3 C 3 . (61) 



The redshift z can be defined by 
which simplifies to 



l + z = ^, (62) 
a 



*-*o = ~7 ^TT + — - 1 )*o + -7^ yrT 1 + z ~- 63 

w{e — 3)Ho V w / w(3 — e)H 

Since, w > and e < 3 always, in our assumption. So for early universe i.e., for z — » oo we 
get 

*-*«>« , * +(-~l)to, (64) 
and for late universe i.e., for z — > — 1 we get 

2" , , ™f§=«2 , x 

*- t -so»^ (1+1,) " T - (65) 

B. Proper distance 

If a photon emitted by a source and received by an observer at time t then the proper 
distance between them is defined by 

[ a °da f to dt 
d = a — = a — , (66) 
Ja aa J t a 

which simplifies to 

, 2-(l + ^^ 2 

(2 + w(e-3))H (2 + w (e-3))H - 1 ' 

C. Luminosity distance 

If L be the total energy emitted by the source per unit time and £ be the apparent 
luminosity of the object then the luminosity distance is defined by 

d *=(ib)'- < 68 > 

Now for our model, 

u, \ 2 W0 (l + z) 2+ ^ 2(1 + z) 
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D. Angular diameter 

The angular diameter of a light source of proper distance D observed at t is defined by 

S-W + 'f. (70) 
d L 

The angular diameter distance (cIa) is defined as the ratio of the source diameter to its 



angular diameter as 



which is simplifies to 



d A =j = d L (l + zy 2 , (71) 



d _ 2-° + 2(l + z)-i 

A (2 + w(e-3))H (2 + w (e-3))# ' 1 1 



The angular diameter distance is maximum at 



*- = i . 2+ : e r 3 ' ) - 1. (73) 



iu(3-e) 2 + w (e-3) 
and hence the maximum angular diameter distance will be 



w(e-3) 



Ml 2 J-%>\~ ( 2 + ro (e-3) 

(<J*)«, = ^ (2 J U(3 -0(2 + ^-3)) J ■ (?4) 



V. CONCLUSIONS 



In this work, we have considered a cosmological model of the homogeneous and 
isotropic FRW universe with variable Newton's gravitational constant and cosmological 
constant. The solutions have been obtained for flat model with barotropic fluid and some 
particular form of cosmological constant (i.e., A = eH 2 ). The cosmological parameters and 
deceleration parameter have been obtained for dust, radiation and stiff perfect fluid. The 
statefinder parameters are analyzed for these three types of fluids and have shown that 
these depends only on w and e. From figure 1, we see that q decreases with w decreases 
and e increases. Also figures 1-3 describe the natures of statefinder parameters for these 
three types of fluids. Further the lookback time, proper distance, luminosity distance and 
angular diameter distance have also been investigated. The maximum angular diameter 
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distance has also been found in our model. 



[1] P. S. Wesson (1978), Cosmology and Geophysics (Oxford : Oxford University Press ); P. S. 

Wesson (1980), Gravity, Particles and Astrophysics ( Dordrecht : Rieded ). 
[2] R. K. Tiwari, Astrophys. Space Sci. 321 147 (2009); M. Jamil et al, Eur. Phys. J. C 60 149 

(2009); M. Jamil et al, Phys. Lett. B 679 172 (2009); M. R. Setare and M. Jamil, J CAP 1002 

010 (2010). 

[3] P. A. M. Dirac, Proc. R. Soc. A 165 119 (1938) ; 365 19 (1979) ; 333 403 (1973); The General 

Theory of Relativity (New York : Wiley) 1975. 
[4] F. Hoyle and J. V. Narlikar, Proc. R. Soc. A 282 191 (1964); Nature 233 41 (1971); C. Brans 

and R. H. Dicke, Phys. Rev. 124 925 (1961). 
[5] Y. B. Zeldovich, Sov. Phys.- JETP 14 1143 (1968); Usp. Fiz. Nauk 11 384 (1968); P. J. E. 

Peebles and B. Ratra, Astrophys. J. 325 L17 (1988). 
[6] P. G. Bergmann, Int. J. Theor. Phys. 1 25 (1968); R. V. Agoner, Phys. Rev. D 1 3209 (1970). 
[7] A. Banerjee, S. B. Dutta Chaudhuri and N. Banerjee, Phys. Rev. D 32 3096 (1985); O. 

Bertolami Nuovo Cimento 93B 36 (1986); Fortschr. Phys. 34 829 (1986); Abdussattar and 

R. G. Vishwakarma, Class. Quantum Grav. 14 945 (1997). 
[8] D. Kaligas, P. Wesson and C. W. F. Everitt, Gen. Rel. Grav. 24 351 (1992). 
[9] A - M. M. Abdel Rahaman, Gen. Rel. Grav. 22 655 (1990); M. S. Berman, Gen. Rel. Grav. 

23 465 (1991); A. Beesham, Int. J. Theor. Phys. 25 1295 (1986). 
[10] Weinberg, S. (1972). Gravitation and Cosmology, (Wiley, New York). 
[11] Abers, E. S., and Lee, B. W. (1973). Phys. Rep. 9, 1. 
[12] Langacker, P. (1981). Phys. Rep. 72, 185. 

[13] A. S. Al-Rawaf and M. O. Taha, Gen. Rel. Grav. 28 935 (1996). 

[14] Arbab I. Arbab, Class. Quantum Grav. 20 93 (2003); Astrophys. Space Sci. 291 141 (2004). 

[15] J. C. Carvalho, J. A. S. Lima and L. Waga, Phys. Rev. D 46 2404 (1992). 

[16] K.D. Krori et al, Gen. Relativ. Grav. 32 (2000) 1439. 

[17] V. Sahni et al, JETP Lett. 77 (2003) 201. 

[18] U. Debnath, Class. Quant. Grav. 25 (2008) 205019. 



12 



[19] X. Zhang, Int. J. Mod. Phys. D 14 (2005) 1597. 

[20] H. Wei and R.G. Cai, Phys. Lett. B 655 (2007) 1. 

[21] X. Zhang, Phys. Lett. B 611 (2005) 1. 

[22] J.Z. Huang et al, Astrophys. Space Sci. 315 (2008) 175. 

[23] W. Zhao, Int. J. Mod. Phys. D 17 (2008) 1245. 

[24] M. Hu and X.H. Meng, Phys. Lett. B 635 (2006) 186. 

[25] W. Zimdahl and D. Pavon, Gen. Rel. Grav. 36 (2004) 1483. 

[26] Y. Shao and Y. Gui, Mod. Phys. Lett. A 23 (2008) 65. 

[27] W. Chakraborty and U. Debnath, Mod. Phys. Lett. A 22 (2007) 1805. 



[28] S. Li et al, |arXiv: 100 2. 3867 [astro-ph.CO] 

[29] T. Singh et al, Gen. Relativ. Grav. 30 (1998) 573. 



[30] A. I. Arbab, astro-ph/9810239 



13 




FIG. 1: The deceleration parameter q is plotted for different choices of e and w. 
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FIG. 2: The statefinder parameters are plotted for w = 1 and e = 1...3. 




FIG. 4: The statefinder parameters are plotted for w = 2 and e = 1...3. 



